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Abstract. Effects of light millicharged dark matter particles on primordial nucleosynthesis 
are considered. It is shown that if the mass of such particles is much smaller than the electron 
mass, they lead to strong overproduction of Helium-4. An agreement with observations can be 
achieved by non- vanishing lepton asymmetry. Baryon-to-photon ratio at BBN and neutrino- 
to-photon ratio both at BBN and at recombination are noticeably different as compared to 
the standard cosmological model. The latter ratio and possible lepton asymmetry could be 
checked by Planck. For higher mass of new particles the effect is much less pronounced and 
may even have opposite sign. 
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1 Introduction 

During the last year there arose renewed interest to the impact of possible light dark matter 
particles on Big Bang Nucleosynthesis (BBN) and Cosmic Microwave Background radiation 
(CMB) [1-3]. The interest was stimulated by an observation that the number of the effective 
neutrino species at BBN is possibly larger than that predicted by the standard model, N e ff = 
3.046 (see e.g. Ref. [4]). Indeed, according to the analysis of the Helium-4 abundance 
N e ff = 3. 7lo*y [5]- This unknown relativistic form of matter got the name dark radiation. 

The effect of additional particles on light element abundances is model dependent and 
may lead both to a rise and to a decrease of N e ff. This has been already found in earlier 
works with MeV dark matter particles [6, 7]. It was argued [7] that if the new particles, X, 
are in thermal equilibrium with neutrinos they would lead to an increase of the Helium-4 
fraction for mx < 10 MeV, while if they are in equilibrium with electrons, positrons and 
photons their effect is opposite. In the first case they correspond to an increase of N e ff and 
in the second case N e ff becomes lower. 

In this paper we also consider an impact of new light particles on BBN. It is assumed 
that A-particles have a small electric charge, much smaller than the charge of electron, 
e'<e. We consider both cases of mx < m e and mx > rn e . The model with mx > rn e has 
been considered in detail previously [7] . In this paper we concentrate on the case mx < n^e ■ 
As we shall see in what follows, during the BBN epoch these dark matter particles were in 
thermal equilibrium with e + , e~, and photons and decoupled from neutrinos. Annihilation 
of e + e~ pairs proceeds both into photons and into dark matter pairs, XX. Later on the 
bulk of dark matter particles annihilates into photons. Here we consider the model where 
this happens after BBN, but of course this should occur long before hydrogen recombination 
to avoid strong constraints on CMBR spectrum distortions. Light elements abundances in 
this model are altered because of: 

i) Increased Universe expansion rate during BBN; 



- 1 - 



ii) Increased neutrino temperature with respect to the photon temperature because the en- 
tropy released in e + e~ -annihilation is now distributed between photons and X-particles; 

hi) Increased baryon-to-photon ratio during BBN. 

We will treat the influence of these factors on BBN analytically, and then calculate numer- 
ically in Sec. 3. But first let us define the presently allowed range of parameters of the 
model. 

2 Existing constraints on milli-charged particles 

We start with the review of existing bounds on the mass and electromagnetic coupling, a 1 , 
of milli-charged particles, in the range were they can be a (part of) dark matter. 

2.1 Laboratory bounds 

For mx < m e the best particle physics bound on a' follows from the data on positronium 
decay to three photons [9], according to which e' < 3.4- 10 -5 e, that is an = 10 11 (e') 2 /47r < 1, 
see also ref. [10], where the similar limit was obtained. For very light particles, mx < 1 keV, 
the best experimental limit e' < 10 -5 e comes from the reactor experiments, [11]. Heavier 
X-particles, mx > m e , are less restricted, e.g. for mx = 100 MeV the bound is a' < 

3 • 10~ 9 [10]. Stringent constraints on a MeV scale dark matter are provided by the fixed 
target experiments, in particular by the LSND [12], see Ref. [13] for a detailed discussion. 
Corresponding bounds may not be applicable to our model because extra light vector boson 
is not necessarily present and we do not aim to explain 0.511 keV annihilation line. For the 
milli-charged scalar particles the main reaction to consider is 7r° — > 7+JJ. Branching of this 
process will be of order a', i.e. the number of produced X particles will be a' times the number 
of produced neutrino. X particles scatter in the detector with the cross-section ~ act / E 2 . 
It will be comparable to the weak cross-section at (E/GeV) 4 ~ 10 10 a(a') 2 ~ 10~ 14 a 2 1 , i.e. 
at E < MeVy / aii/3, which looks safe. 

2.2 Constraints from stellar evolution 

The consideration of stellar cooling also allows to restrict the interaction strength of light 
X-particles [14]. These particles could be abundantly produced in stellar cores and if their 
mean free path is shorter than the stellar radius, they would efficiently carry out stellar 
energy, strongly amplify the stellar luminosity, and diminish the star age. Inside the Sun 
they could be created in the process ee — > eeXX. The cross-section of this reaction can be 
roughly estimated as a(ee — > eeXX) ~ a 3 c//m 2 . The energy transferred to X-particles in 
the solar core per unit time is about 

[internal „ rfj^E^ ~ l 44 a n erg/s , (2.1) 

where n e = 6 • 10 25 /cm 3 is the density of electrons in the solar core, R c ~ 10 10 cm is the 
core radius, Ex ~ keV is the X energy, and v e ~ (T c /m e ) ~ 0.05 is the thermal velocity of 
electrons in the solar core. However, the mean free path of X-particles is much shorter than 
R c and their emission is suppressed by the factor (lf ree /R c ) 2 . The mean free path due to 
eX-scattering can be estimated as 

Ifree = iPeX^X^e) , (2.2) 
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where a e x = 4:iraa' 'L/ '(m x v x ) * s ^ ne Coulomb-like transport cross-section, and L is the 
Coulomb logarithm. Taking mx = 0.3 keV and vx ~ 1, we find lf ree ~ (4/L) cm. Cor- 
respondingly the solar luminosity due to radiation of X-particles would be below 10~ 8 Lq, 
where Lq = 4 • 10 33 erg/sec is the solar luminosity. 

The X-luminosity of white dwarfs (WD) can be estimated similarly. The electron 
density inside WD is 10 4 larger than the solar one, while the radius is 10 times smaller than 
R c . So L™ ternal would be 10 5 times larger. However the mean free path of X-particles would 
be 10 3 times smaller. The factor 10 -4 comes from the larger n e and factor 10 comes from 
smaller eX-cross-section because of larger X energy. So finally the X-particle luminosity of 
white dwarfs would be about 10 26 erg/sec, which looks safe. 

There is a competing process of X-production for WD, namely decay of plasmons into 
the XX-pairs, because the plasma frequency in WD is about 10 keV and so it may be larger 
than mx- However this process is less efficient than X-production in the above discussed 
reaction or at most comparable. Similar situation takes place in red giants. 

Note that surprisingly strong limit a' < 10 -30 for mx < 10 keV was obtained from 
the stellar cooling in ref. [15]. However, this bound is valid if a' < 10~ 18 is granted. The 
authors of ref. [15] assumed the validity of the last bound using BBN considerations. In our 
approach this BBN bound is invalid since we look if it is possible to allow millicharged dark 
matter particles at the expense of some additional modification of BBN, e.g. by introduction 
of neutrino asymmetry. 

2.3 Cosmological bounds 

Late annihilation of X-particles, when they go out of chemical equilibrium with photons may 
distort the energy spectrum of CMB. The noticeable deviations from equilibrium occurred at 
T ~ mx/W, when the bulk of X-particles was annihilating. At this temperatures the energy 
density of X-particles is of the order of 4 • 10 -2 Vtx / m\ of the energy density of the CMB pho- 
tons. Kinetic equilibrium of energetic photons, created in XX-annihilation, is restored by the 
elastic 7 e-scattering very quickly, with the characteristic time of approximately 4(keV/T) 3 
sec. However, chemical equilibration of photons demands inelastic reaction 7 e — > 27 e whose 
probability is approximately five orders of magnitude smaller, so the effective reaction time is 
about 4 • 10 5 sec at T = 1 keV. It is quite close to the cosmological time at this temperature. 
As is known [16], large non-equilibrium energy influx, of the same order of magnitude as 
the energy density of CMB, into the cosmological plasma would be perfectly thermalized if 
it took place at z > 10 . Smaller influx could be thermalized at lower T, so the chemical 
potential of the CMB photons could be presumably reduced below the observational bound 
fi < 10~ 4 T. 

The fraction of X-particles which were annihilating during recombination at T ~ 0.2 eV 
equals to 0.2 eV/(m x /10) ~ 10~ 3 (for S-wave and even smaller for P-wave annihilation). 
Therefore, the energy fraction of keV photons will be 10 -3 times the ratio of energy densities 
of X-particles and photons during this epoch. If Qx is> sa y ; 0.1 of the baryonic contribution 
to £1, then the corresponding perturbation of the chemical potential will not exceed 10 -4 , and 
therefore will not exceed observational bounds. Taking into account that spectrum distortions 
appear only after photon degradation in energy, down from keV to the CMB energy, we see 
that the fraction of X-particles can be even higher. 

Another very interesting and important bounds are coming from the analysis of CMBR 
anisotropies in the presence of millicharged particles [17]. One concludes that the fraction 
of such particles cannot be too large at the recombination if they are coupled to radiation. 
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This bound is not really restrictive in our situation since we are not assuming here that 
millicharged particles constitute significant fraction of matter presently. 



3 BBN in the presence of milli-charged particles 

3.1 Very light dark matter, rax < m e 
3.1.1 Cosmological abundance 

Light charged particles will be in thermal equilibrium with electromagnetic plasma in the 
early Universe when 

a{e + e~ -> XX) n e > H, (3.1) 

where 

V 90 m pi m pi y ' 

is the Hubble parameter at the radiation dominated cosmological stage and n e ~ T 3 is the 
number density of electrons. Parameter = 10.75 + (7 / '4){N e f f — 3) counts the number of 
particle species at m e < T < m^. The cross-section of XX production by electron-positron 
pairs for high electron energy, E e > m e , is 

a(e + e- XX) ~ aa'/T 2 , (3.3) 

where a 1 = (e') 2 /47r. So the process is in equilibrium roughly at m e < T < aa'm pi and we 
expect that at T ~ 1 MeV X-particles have the equilibrium energy density, corresponding 
to one neutrino species up to the Bose- Fermi factors. 

The late time (frozen) cosmological density of X-particles is determined by their anni- 
hilation into two photons and according to the standard calculations, see e.g. [8], is equal 
to 

c 10n 7 ln[o-(XX 2^)v x m x mpi} 

n x ~ T^vv — T^T^ ' v J - 4 J 

a {XX — > Z^jvxrnx m Pi 

where vx is the center-of-mass velocity of X-particles and 

va(XX -> 2 7 ) = ir(a'f/m x . (3.5) 



So the mass density of X-particles at the present time, if they are stable, should be equal to 

p c x = 10keV/cm 3 ( —\ [1 + 0.12 lnfafi/mi)], (3.6) 
V«ii/ 

where mi = mx/keV and an = a 1 /10 -11 . For comparison, the present day total cosmolog- 
ical energy density is about 5keV/cm 3 . 

Estimating p c x we have neglected the plasmon decay j p i — > XX. This is justified 
because at high temperatures, T > m e , when the plasma frequency, uj p i ~ 0.1T, exceeds mx 
an account of the plasmon decay slightly shifts equilibrium condition for X-particles which are 
in equilibrium anyhow. At smaller temperatures, T < m e , when XX-annihilation is frozen, 
the plasma frequency is either Bolzmann suppressed, ~ exp(-m e /mj), or suppressed by the 
smallness of a' and the plasmon decay may be neglected. 
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3.1.2 Qualitative discussion of BBN 

Let us return now to evolution of X, , and photons at BBN. We assume that AT-particles 
do not have anomalously strong interactions with neutrinos (though this possibility may be 
interesting) and thus at T below 1 MeV neutrinos are decoupled from X and the electro- 
magnetic part of the primeval plasma. Due to reaction e + e~ — > XX and elastic scattering 
X-particles remain in good thermal contact with electrons and positrons. Indeed, the ratio 
of the reaction rate to the Hubble parameter is 

1_ nx_ _ Traa'nl m P i 
H nx 5m^nxT 2 

and it is much larger than unity even for T <C m e , when the electron number density is 
exponentially suppressed. Correspondingly e + e~-annihilation equally well heats up photons 
and A-particles. So the entropy of electron-positron pairs which in the usual case was totally 
transferred to photons, is now distributed between 7, X, and X. Correspondingly the entropy 
factor, which in the standard model was 11/4, now becomes 15/8 (assuming that X have spin 
zero) . Hence the relative energy density of one neutrino species with respect to photons after 
e + e- -annihilation instead of (7/8)(4/ll) 4 / 3 = 0.227 would become (7/8)(8/15) 4 / 3 = 0.378. 
However, this change takes place asymptotically at T <C m e but at T ~ m e the effect is much 
smaller, which can be easily estimated analytically. 

This rise of relative neutrino density has two- fold effect on BBN. First, rising density of 
v e would shift the temperature of neutron-proton freezing to smaller values and thus leads 
to a decrease of N e tt. On the other hand, the rise of the energy density of relativistic 
species increases N e ff through the corresponding increase of the Hubble parameter at the 
same temperature. At the moment of n/p- freezing the first effect is stronger but later on the 
second effect dominates. It is especially pronounced in its impact on the time of the onset of 
BBN. The formation of first light elements started roughly at 

Tbbn « E h / ln(l/7/) « 70 keV, (3.8) 

where n is the ratio of baryon to photon number densities, n = ns/n^. The time when 
Tb bn is reached is determined by the number of relativistic species in the plasma according 
to tBBN ~ l/G/5* Now tsBN becomes considerably shorter and less neutrons have 
time to decay prior to BBN. It can be described as quite large rise of N e ff. Below we 
present numerical calculations of primordial nuclei production for very light A-particles with 
mx < Tbbn which are in very good agreement with our simple analytical estimates. 

The baryon-to-photon ratio at BBN in this model is different from that determined 
from CMBR. The latter, according to WMAP observations, is 77 = 6.19 x 10~ 10 [18]. We 
can determine 770 at the moment of n/p- freezing, which takes place at T ~ 0.7 MeV using 
entropy conservation. Annihilation of N extra degrees of freedom into photons during the 
epoch between BBN and recombination dilutes 77 by the factor of (2 + N)/2. Therefore, 
the baryon-to-photon ratio right after e + e~ annihilation should be larger by this factor as 
compared to the standard cosmological model. Annihilation of e + e~ pairs leads to dilution 
factor of (ll+2N)/(4+2N). The product of these two factors determines the initial condition 
for the baryon-to-photon ratio: 

11 + 2N . . 

770 = n. (3.9) 

In particular, with one charged extra scalar N = 2 and 7/0 = (15/4) 77, while in the standard 
model this ratio is = (11/4) rj. One should keep in mind however, that in the standard 
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Figure 1. Left panel: Helium abundance vs lepton asymmetry in our model is shown by solid curve. 
Gray band indicates observationally allowed range of helium abundance. Range of x, required to 
match observations, is shown by vertical dotted lines. Right panel: Deuterium abundance vs lepton 
asymmetry. Gray band indicates observationally allowed range of D/H. Vertical dotted lines are from 
the left panel. 



model r] changes to its asymptotic value determined by CMB in the interval from n /p- freezing 
down to almost complete e + e _ -annihilation, which takes place before T = Tbbn, while in 
our scenario the ultimate entropy release from XX — > 2j takes place at T ~ rax <C Tbbn- 
The ratio of the neutrino temperature to the photon one at n/p-freezing in our model 
remains the same as in the standard case, T^/T 7 = (4/11) 1 / 3 = 0.71, while after e+e~- 
annihilation but prior to Tgpjjsj it rises cis 



2N N 1/3 
11 + 2N 



T - = ( Tm^-r ) T 7, (3-10) 



so for scalar X-particles it would be T^/T 7 = (8/15) 1 / 3 = 0.81. Subsequently after XX- 
annihilation to photons this ratio drops down by [2/(2 + iV)] 1 / 3 and turns to T v /T^ = [4/(11 + 
27V)] 1 / 3 which for N = 2 becomes T u /T^ = 0.64. In other words, the energy density of 
neutrinos determined by CMB would be 1.5 smaller than in the standard model. 

3.1.3 Numerical results 

We have calculated abundances numerically. To this end we have modified accordingly the 
publicly available numerical code of Kawano [19], so that all effects mentioned above are 
automatically taken into account. 

Following Refs. [20, 21] we use the conservative constraints on light elements abundances 

0.240 < Y p < 0.258, (3.11) 

1.2 • 10~ 5 < D/H < 5.3 • 10~ 5 , (3.12) 

0.57 < 3 He/D < 1.52, (3.13) 

7 Li/H > 0.85 x 10~ 10 . (3.14) 

Without lepton asymmetry the He abundance turns out to be in clear conflict with 
observations. In the standard model it is BBN which gives the strongest constraints on 
asymmetry between neutrinos and antineutrinos. Due to strong mixing between different 
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Figure 2. Left panel: 3 He/D ratio vs lepton asymmetry. Gray band indicates observationally allowed 
range of this ratio. Vertical dotted lines are from Fig. 1. Right panel: Li/H ratio vs lepton asymmetry. 
Observationally allowed range is above horizontal dotted line. 



neutrino species chemical potentials in all neutrino sectors are the same and in the standard 
model they are bounded from above as x = f^/T < 0.07 [22]. 

In the present case, on the contrary, we can introduce lepton asymmetry to bring Helium 
abundance in accord with observations. 

Calculated He abundance as a function of x is shown in Fig. 1, left panel, for the case of 
charged scalar, N = 2. Colored band represents the range given by inequality (3.11), which 
in turn determines the required range of x shown by vertical dotted lines. 

Corresponding results for D/H, 3 He/D and Li/H ratios are shown in Figs. 1 and 2. As 
we see, ther calculated ratios are marginally consistent with conservative constraints (3)-(6). 
We should stress that 7 Li abundance in our model is about 50% above the standard BBN 
predictions. This worsens the standard model problem with 7 Li and we need even more 
stellar depletion that is typically invoked. 

3.2 Heavier dark matter, mx > m e 

Situation becomes very much different for heavier X-particles, especially if mx > rn e , and was 
studied in detail in Ref. [7]. In this case the cosmological number density of X-particles would 
be determined by their annihilation into e + e~ . Taking for the annihilation cross-section the 
approximate expression 

va(XX — > e + e~) = koloit /m\, (3.15) 
where r = yj\ — Am1/m 2 x and assuming for simplicity that mx 3> m e , we find: 

p\ « 2 • 10- 2 ^ (^) [1 + 0.05 ln( ail /m 3 )}, (3.16) 
cm \atu J 

where 771,3 = m^/MeV. 

If 777-x is sufficiently small, such that the equilibrium number density of X-particles is 
non-negligible at n/p- freezing, then their presence would speed up the cosmological expansion 
and in this sense is equivalent to some dark radiation, though X-particles at this stage were 
non-relativistic or at most semi-relativistic. On the other hand, additional e + e _ -pairs from 
XX-annihilation would diminish relative contribution of neutrinos with respect to photons 
and this results in a decrease of N e tt. For example for mx — m e the ratio of neutrino to 
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photon temperature would be (T^/T 7 ) 3 = 4/(11 + 2N) and if N = 2 we obtain T u /T y = 0.64 
already at BBN but not much later as it was in the case of light X-particles. So depending 
upon rax the overall effect of X-particles on BBN may be of either sign. E.g., as it was 
shown in Ref. [7], the existence of millicharged particles in the mass range mx = 4-10 
MeV can even improve the overall agreement between the predicted and observed 2 H and 
4 He abundances. 

4 Conclusion 

We have found that assuming non-zero lepton asymmetry one may avoid the BBN upper 
bounds on the charge of possible light milli-charged particles. This in turn opens a way to 
modify the standard BBN predictions and to mimic possibly observed dark radiation. An 
interesting and testable feature of the model is that the values of the baryon and lepton 
asymmetries (i.e. the ratios n^/n 7 and n^/n 7 ) are different at BBN and recombination 
epochs. 

The model discussed here was stimulated by our work on the role which X-particles 
might play in the generation of large scale magnetic fields after hydrogen recombination [23] . 
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